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Abstract 

We discuss the relation between cosmic-ray induced interstellar gamma-ray and radio emission from our Galaxy, 
and emphasize the importance of their parallel study. We give an overview of results from Fermi-LAT on the inter- 
stellar gamma-ray emission, and then focus on complementary studies of the radio emission from the Galaxy. 
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1. Introduction 

Cosmic rays (CRs) propagate in the Galaxy and inter- 
act with the interstellar medium (ISM). From their inter- 
actions with the gas and the interstellar radiation field 
(ISRF), gamma-ray diffuse emission is produced. The 
mechanisms of this diffuse emission are: neutral pion 
decay from hadronic interactions of CR protons and 
heavier nuclei with gas, inverse-Compton (IC) scatter- 
ing from CR electron and positron interactions with the 
ISRF and the cosmic microwaves background (CMB), 
and bremsstrahlung from CR electron interactions with 
gas (see eg. Q). The same CR electrons and positrons 
producing gamma rays also generate diffuse radio emis- 
sion, via synchrotron radiation in Galactic magnetic 
fields (B -field). Interstellar synchrotron emission ex- 
tends from a few MHz to tens of GHz and is one of the 
major components of radio and microwave band radia- 
tion. This emission depends on the CR electron spec- 
trum and distribution in the Galaxy, and on the B -field 
(see eg. l29l ). 

Observations of gamma-ray and radio emission from 
our Galaxy are an important tool for studying CR, the 
ISRF, interstellar gas and the B -field. Using the combi- 
nation of these energy bands and direct measurements 
of local CRs, we can gain information on CR propaga- 
tion models, CR source distribution, and CR spectra in 



interstellar space. The importance of parallel studies of 
radio and gamma-ray emission was recognized already 
in the 80's, eg. in 031, d, E3 and (B) and more re- 
cently in [18], 1 29]. This approach has advantages over 
studying these data separately, since the CR electrons 
are traced by both gamma-ray and radio emissions and 
so degeneracy is removed. CR electrons are measured 
in direct observations, but these are affected by solar 
modulation at lower energies. 

We present here an overview of the latest results with 
Fermi-LAT on the gamma-ray diffuse emission induced 
by CR nuclei and electrons. Then we focus on the 
complementary studies of the synchrotron emission in 
the light of the latest gamma-ray results. These stud- 
ies made use of the numerical model of CR propagation 
and interactions in the Galaxy, GALPROP 031 [MS, de- 
scribed in Section 3. 

2. Cosmic rays throughout the Galaxy 

CRs are believed to originate in supernovae rem- 
nants (SNRs). During propagation in the Galaxy, the 
CR spectra are steepened by energy losses and energy- 
dependent diffusion, but can also be affected at low en- 
ergies by diffusive reacceleration. The spectrum orig- 
inating in the sources is called the injected spectrum, 
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the propagated spectrum pertains to interstellar space, 
while that measured in the Solar System is the solar- 
modulated spectrum. Inverse-Compton scattering, syn- 
chrotron radiation and ionization are responsible for en- 
ergy losses of electrons and positrons. Direct mea- 
surements of CRs by balloon experiments and satellites 
constrain models of propagation and local CR spec- 
tra. However solar modulation complicates the inter- 
pretation of the observations at energies below about 10 
GeV/n. An important point is that gamma-ray and syn- 
chrotron emission probe the interstellar CR spectra free 
from the effects of modulation. 



3. Modelling interstellar emission with GALPROP 

GALPROP is a software package for modelling the 
propagation of CR in the Galaxy and their interactions 
in the ISM. Descriptions of the GALPROP software can 
be found in fT5lL 1261 , l27l and l34l and references 
therein. This project started in the late 1990's |[T4ll24ll . 
Since then, it has been continuously improved and up- 
dated. It enables simultaneous predictions of obser- 
vations of CRs, gamma rays and recently synchrotron 
radiation (Tgl OH (29). GALPROP models were used 
to analyse the diffuse gamma rays detected by CGRO 
EGRET and COMPTEL, INTEGRAL ED, and more 
recently radio surveys and WMAP data [19] |29) . It is 
used in the Fermi-LAT collaboration for interpreting the 
physics of the Galactic emission (6). For details we re- 
fer the reader to the relevant papers |6, 29] and the ded- 
icated websit^] GALPROP solves the transport equa- 
tion for a given CR source distribution and boundary 
conditions, for all required CR species. It takes into ac- 
count diffusion, convection, energy losses, and diffusive 
reacceleration processes. Secondary CRs produced by 
collisions in the ISM and decay of radiative isotopes are 
included. The propagation equation is solved numeri- 
cally on a user-defined spatial grid in 2D or in 3D, and 
energy grid. The solution proceeds until a steady-state 
solution is obtained, starting with the heaviest primary 
nucleus, and then electrons, positrons, and antiprotons 
are computed. GALPROP models gamma-ray emission 
of the three components, pion decay, bremsstrahlung 
and Inverse Compton for a user-defined CR source dis- 
tribution and CR spectra. Gas maps and ISRF are 
provided and updated to the most recent observations 
l6l . More recently, GALPROP calculation of interstel- 
lar synchrotron emission has been improved |29l and it 
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Figure 1: Global luminosity spectra of the Galaxy. Line styles: inter- 
stellar radiation field, including optical and infrared (magenta solid), 
protons (red dotted line), helium (blue dotted line), primary electrons 
(green dotted line); CR-induced diffuse emissions (solid lines): IC 
(green), bremsstrahlung (cyan), n decay (red), synchrotron (black), 
total (blue). CR particle luminosities are also shown as labelled. 



has been extended to include also synchrotron polariza- 
tion, and free-free emission and absorption lfT9l . 

In (28) we used the GALPROP code to derive the 
broad-band CR-induced luminosity spectrum of the 
Galaxy. Figure [T] shows an example of the modelled 
luminosity of the Galaxy. The total luminosity in syn- 
chrotron and inverse Compton emission are similar in 
magnitude. We found the Galaxy to be a fairly good lep- 
ton calorimeter, when both inverse Compton and syn- 
chrotron losses are accounted for. 



4. Interstellar gamma-ray emission 

Interstellar gamma-ray emission accounts for at least 
60% of the gamma-ray photons detected by Fermi- 
LAT. An early result from Fermi-LAT (TJ was a study 
of this gamma-ray emission. This study was car- 
ried out using data from intermediate Galactic lati- 
tudes, where the interstellar emission comes mainly 
from local CR interactions with the ISM. The spec- 
trum was found to be well described by pion decay, 
IC and bremsstrahlung, plus the contribution from point 
sources and an isotropic component, which included the 
extragalactic background. This result was confirmed in 
l2ll and Q, where the local HI gamma-ray emissivity 
was derived for the 2nd and 3rd Galactic quadrants. The 
latter studies also showed an increase in the H 2 -to-CO 
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factor (Xco) with Galactocentric radius and larger in- 
tensities in the outer Galaxy with respect to the a priori 
model based on local CR measurements. A flatter dis- 
tribution of CR sources in the outer Galaxy or a larger 
amount of gas could explain this effect. Models used for 
these papers, generated with GALPROP, included diffu- 
sive reacceleration with a 4 kpc halo height, the locally- 
measured CR proton and Helium spectra, and the elec- 
tron spectrum as measured by Fermi [[3) and (4)1 . Models 
used Alfven velocity in the range 30-39 km s" 1 (i.e. see 
ED, ED, ESI, E3, (H). CR sources were assumed 
to follow supernova remnants as traced by pulsars. A 
preliminary study including the Galactic plane was pre- 
sented in [32], where the spectrum of the inner Galaxy, 
longitude and latitude profiles were shown. The model 
reproduced the Fermi data over the sky within 20%; the 
inner Galaxy revealed an excess at GeV energies that 
can plausibly be explained by unresolved sources. It 
was shown that latitude profiles improved with a larger 
halo height of about 10 kpc, consistent with the pre- 
vious works mentioned above. More recently (6) in- 
vestigated models of the gamma-ray diffuse emission, 
taking into account uncertainties associated with the as- 
trophysical parameters. The distribution of cosmic-ray 
sources, the size of the cosmic-ray halo, and the inter- 
stellar gas distribution were varied within observational 
limits. They produced a grid of GALPROP models and 
compared them with Fermi-LAT data for the first 21 
months of observations. They found a general agree- 
ment between models and Fermi-LAT data, showing 
that the physics of the diffuse emission is broadly un- 
derstood. The gamma-ray fits improved using dust as a 
tracer of gas, confirming the presence of molecular hy- 
drogen that is not traced by CO: the so-called dark gas 
as found by ifTOl and recently confirmed by [20]. No 
single model gave a best fit over all sky regions. In gen- 
eral, the models were consistent with the data at high 
and intermediate latitudes. However, all models under- 
predicted the data in the inner Galaxy for energies above 
a few GeV. In the outer Galaxy, they found that the data 
prefer models with a flatter distribution of cosmic-ray 
sources or greater gas density. In general, all models 
were within 15% of data, except for some regions such 
as Loop I, and the structures coincident with the bub- 
bles found by [30], and the above-mentioned excess in 
the outer Galaxy. 

5. Interstellar radio emission 

Interstellar synchrotron radiation is the major contrib- 
utor to the Galactic radio emission from tens of MHz 



to a few GHz. At higher frequencies and in the mi- 
crowave bands, free-free and dust emission tend to dom- 
inate. Synchrotron emission is also a foreground com- 
ponent for CMB studies. Recent studies of the Galactic 
synchrotron emission using CR electron measurements 
were published in (I7l[29]|. Below a few GeV the local 
interstellar electron spectrum cannot be directly mea- 
sured, because CR electrons with energy lower than a 
few GeV are affected by solar modulation. The syn- 
chrotron spectral index depends on the CR electron and 
positron spectral index, while synchrotron intensity de- 
pends on the B-field intensity and electron flux. In par- 
ticular CR electrons from 500 MeV to 20 GeV produce 
synchrotron emission from tens of MHz to hundreds of 
GHz for a B-field of few //G, and hence this can be 
used in conjunction with direct measurements to con- 
struct the full interstellar electron spectrum from GeV 
to TeV 

Using a collection of radio surveys and WMAP data, 
out of the Galactic plane, |29l found that the propagated 
local interstellar electron spectrum turns over below 
a few GeV, with an ambient interstellar index around 
2, independent of modulation. We tested propagation 
models, generated with GALPROP, in order to constrain 
the injected CR electron spectrum, before propagation. 
We found the need for an injection electron spectral in- 
dex between 1.3 and 1.6 below 4 GeV (see Fig[2]). One 
result was the recognition of the importance of includ- 
ing secondary positrons and electrons in synchrotron 
models. While plain diffusion models fitted the data 
well, standard reacceleration models often used to ex- 
plain CR secondary-to-primary ratios were not consis- 
tent with the observed synchrotron spectrum, since the 
total intensity from primary and secondary leptons ex- 
ceeded the measured synchrotron emission at low fre- 
quencies, as shown in Fig [2] We confirm here that plain 
diffusion models or diffusive reacceleration models with 
small Alfven velocity are preferred to existing reaccel- 
eration models, as found previously in [29]. Spectra 
of the total synchrotron emission for high latitudes re- 
gions and for the case of lower Alfven velocity than 
usually assumed are shown in Fig [2] Reacceleration 
models with Alfven velocity around 20 km s" 1 also fit 
the synchrotron data. A detailed description of our re- 
cent improvements in the radio modelling is given in 
fT9l . An example of our modelling of the radio emis- 
sion is shown on Figure [3] Here spectra of total (I) and 
polarized (P) emission from radio surveys and WMAP 
data are compared with our model in the inner Galaxy, 
showing reasonable agreement. Polarized synchrotron 
emission allows the regular component of the B-field to 
be studied, complementing the total synchrotron which 
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probes the sum of random and regular fields. Models 
of the regular B -field from the literature can be investi- 
gated, using CR propagation models based on recent CR 
measurements. We studied the sensitivity of the syn- 
chrotron model to various formulations of the regular 
B -field in the literature (see caption to Fig [4] for details 
of the models), finding that while all models resemble 
the data to first order, no single model of the ones used 
fitted best all the data for the whole sky. Skymaps at 
408 MHz and polarized emission at 23 GHz are shown 
in Fig|4]for various models of B-field, halo size, and CR 
source distribution. Note that local features like Loop I 
are not included in the model, and these account for 
much of the additional structure seen in the data. The 
models used reproduce the peak in the direction of the 
inner Galaxy for both total and polarized radio compo- 
nents. For the magnetic field models used, 3-4 yuG and 
5-6 yuG respectively for the regular and random compo- 
nent of the local B-field described the data best. This 
regular field is a factor 1.5-2 larger than in the original 
models. We find that increasing the height of the halo 
from 4 kpc to 10 kpc gives better agreement with radio 
data for high latitudes. However it somewhat overes- 
timates the emission at intermediate latitudes. Hence 
no conclusions can be drawn at this point, and the halo 
height should be investigated in more detail in the fu- 
ture. The effect of varying the CR source distribution 
has also been investigated, indicating that a flat distri- 
bution in the outer Galaxy reproduces better the data. 

We model also the expected synchrotron spectral in- 
dex and its variation in the Galaxy. Figure [5] shows an 
example of a skymap of the spectral index from 408 
MHz to 23 GHz, for a given model of magnetic field. 
We see that the variations in the Galaxy are very small 
compared for example with the one obtained by [ 7 ] with 
WMAP where the spectral index range is 2.3-3.4. How- 
ever the trend is similar, with a smaller spectral index 
in the plane, and increasing away from the plane. This 
means that the physical model of electron propagation 
is supported. 

6. Summary 

We have described some recent results on interstellar 
emission related to Galactic CR. We have also presented 
the recent improvements on our modelling of the inter- 
stellar radio emission. The aim of this paper is to high- 
light what we can learn from a parallel study of gamma 
rays and radio emission in a multi- wavelength approach. 

Generally the physics of these interstellar emissions 
is well understood. Models with a CR source distribu- 
tion rather flat in Galactocentric radius beyond 10 kpc 



Figure 2: Synchrotron spectra (blue lines) with primary low-energy 
electron injection index 1.6 for pure diffusion model (upper), for dif- 
fusive acceleration model as in fp (center) and for diffusive reacceler- 
ation model with Alfven velocity 20 km s _1 (lower) at high latitudes. 
Data (red triangles) are from radio surveys and WMAP. 
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Figure 3: Radio spectra of total intensity I (upper) and polarized in- 
tensity P (lower) for the inner Galaxy. B -field intensities are scaled 
with respect to the original models to agree with data. The plots show 
various model components: synchrotron (blue), thermal dust + spin- 
ning dust (magenta upper plot), thermal dust (magenta lower plot), 
free-free (green), free-free+synchrotron (red line) and total (cyan). 
Data are from radio surveys (black triangles) and WMAP (red tri- 
angles). 

reproduce better both gamma-ray and radio data in the 
outer Galaxy. A halo size larger than 4 kpc is indicated 
in both energy bands in order to fit the high latitude re- 
gions. However, no single model is found to reproduce 
the data best over the sky. Present diffusive reacceler- 
ation models, with parameters often used to reproduce 
secondary/primary ratios, are not consistent with syn- 
chrotron data. In order to reproduce synchrotron data, 
plain diffusion models or models with small diffusive 




Figure 4: Brightness temperature skymaps at 408 MHz (left) and po- 
larized 23 GHz (right). Top to bottom are the following models: - 
ASS+RING B-field model based on (31], for CR source distribution 
from 1 28 1, for 4 kpc halo size; - ASS B-field model based on 1 22], for 
CR source distribution from (281 , for 4 kpc halo size; - BSS B-field 
model based on f22l , for CR source distribution from (28], for 4 kpc 
halo size; ASS+RING B-field model based on |31|, CR source distri- 
bution from (28], for 10 kpc halo size; - ASS+RING B-field model 
based on (31], CR source distribution from fl3l . for 4 kpc halo size. 
The intensities of the local magnetic field are increased of a factor of 2 
with respect the original models in order to fit the radio data (see text 
for more details). For comparison for all the models shown the in- 
tensity of the local regular component is 4 //G and for the random is 6 
/jG. The skymaps at the bottom are the data from [TT] and WMAP 1 8 1. 
All skymaps (for each frequency) have the same scale, with Galactic 
longitude 1 = at the center. Units are K (left) and mk (right). 

reacceleration are needed. Gamma-ray studies alone 
leave unanswered questions which can be addressed by 
investigating the radio band, hence putting more con- 
straints on the models. For example, gamma-ray data 
are not very sensitive to the electron spectrum, because 
gamma emission is dominated by the hadronic compo- 
nent rather than the IC emission. In this case, gamma 
rays and radio data give complementary information. 

These studies are also of interest for upcoming data 
from the Planck mission. 
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Figure 5: An example of the modelled synchrotron spectral index 
from 408 MHz to 23 GHz for the full sky. The map is in Galactic 
coordinates. 
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